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The overlap operator is just the simplest of a class of Dirac operators with an 
exact chiral symmetry. I demonstrate how a general class of chiral Dirac operators 
can be constructed, show that they have no fermion doublers and that they are 
all exponentially local, and test my conclusions numerically for a few examples. 
However, since these operators are more expensive than the overlap operator, it is 
unlikely that they will be useful in practical simulations. 
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1 Introduction 

For many years, it seemed that simulating chiral symmetry on the lattice was 
impossible, because of the Nielson-Ninomiya theorem pQ, which, in one form, 
states that it is impossible to have a Dirac operator which is local, has the cor- 
rect continuum limit, hyper-cubic symmetry, no doublers, and anti-commutes 
with 75. A method of avoiding the Nielson-Ninomiya theorem was suggested 
shortly afterwards, in 1982, by Ginsparg and Wilson [2], who proposed that 
the smoothest way of breaking chiral symmetry on the lattice in the contin- 
uum limit is to use a Dirac operator which, rather than anti-commuting with 
75, satisfies the Ginsparg- Wilson relation, 



I will call R the Ginsparg- Wilson function, and it plays a key role in this 
paper. 
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The Ginsparg- Wilson relation was soon forgotten as no solutions were known. 
After this relation was rediscovered in 1998 [3], Martin Liischer showed that 
it implied an exact symmetry of the QCD Lagrangian for any lattice Dirac 
operator obeying the Ginsparg- Wilson relation with a local R, and that this 
symmetry reduces to chiral symmetry in the continuum limit [I]. Inspired by 
the work of Kaplan [5] and using an infinite number of fermion fields to avoid 
the Nielson-Ninomiya theorem, Neuberger and Narayanan had already, a few 
years earlier, found a Dirac operator which satisfied the Ginsparg- Wilson re- 
lation with R — 1, namely the overlap operator [BJ, which can be written in 
the form given in section [2], equations (fT3"j) and f[T^|) . Non-local solutions to 
the Ginsparg- Wilson equation have also been found, for example in [8], and 
the various fixed point actions, including the classically perfect action, [9] are 
known to satisfy the Ginsparg Wilson relation [3] , but in practice a truncated 
form of the Dirac operator has to be used, meaning that the chiral symmetry 
becomes inexact. Fujikawa has constructed the algebraic solutions to a par- 
ticular form of a generalised Ginsparg- Wilson equation [11]. The construction 
of these operators requires two steps: firstly constructing an intermediate op- 
erator from the matrix sign function of a Wilson-type operator, then taking 
a root of that operator. Thus, these operators effectively require a nested se- 
ries of roots of a matrix, meaning that they will be considerably slower than 
overlap fermions in practice. Werner Kerler has constructed generalised chiral 
lattice Dirac operators by solving a particular form of the Ginsparg Wilson 
relation [13], and I discuss the relationship between his work and my own 
in section The overlap operator, Fujikawa's and Kerler's solutions are all 
particular forms of my most general solution. I shall not discuss the generali- 
sation of the Ginsparg- Wilson equation obtained by using different kernels in 
the matrix sign function. 

Of course, not every possible Ginsparg- Wilson function R will lead to a lattice 
chiral symmetry. An easy, though over simplified, way of seeing this is to note 
that if, in an expansion in the lattice spacing, R is inversely proportional 
to the lattice spacing, a, (or worse) then the right hand side of equation 
(CQ) will remain constant (or diverge) rather than reduce smoothly to zero in 
the continuum limit. More strictly, the lattice chiral symmetry requires that 
{75, R} is non- vanishing and local [IB]. This places a number of restrictions on 
the lattice Dirac operator D, such as that if it satisfies the Ginsparg- Wilson 
symmetry and has no doublers then it cannot be ultra-local [17], but at best, 
like the overlap operator [18], exponentially local. 

I define ®c as the class of possible suitable Dirac operators (local, 75-Hermitian, 
with the correct continuum limit, and no doublers) which have their eigenval- 
ues lying on a curve in the complex plane. It is interesting that the overlap 
operator, the generalisations by Fujikawa and Kerler, and the Dirac operator 
in the continuum are all within ®c- in the continuum, the eigenvalue spec- 
trum is a straight line along the imaginary axis; for overlap fermions, the 
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curve is a circle. However, for example, Wilson fermions, which have no chiral 
symmetry, and staggered fermions, which have doublers, are outside Dc- It is 
an interesting hypothesis whether all suitable lattice Dirac operators within 
D c satisfy Ginsparg- Wilson chiral symmetry. In this paper, I demonstrate the 
plausibility of this hypothesis by showing that a class of Dirac operators with 
the correct continuum limit, no fermion doublers, and with their eigenvalue 
spectrum lying on an arbitrary curve in the complex plane are all suitable 
chiral fermions — in that they obey a Ginsparg- Wilson relation with local R 
and are local themselves. 

It is also worth considering the reverse hypothesis, whether all possible Ginsparg- 
Wilson lattice Dirac operators are within the class Dc- In section [521 1 demon- 
strate that this is true for Dirac operators satisfying [D, D*] = 0, a condition 
obeyed in the continuum. However, in general it does not seem to be the case. 
Certain fixed point lattice Dirac operators [31TT9] offer one known counterex- 
ample, and I discuss others in section 12721 However, these Dirac operators are 
related to a Dirac operator within Dc by a simple chirally invariant transfor- 
mation. However, there are distinct homotopy classes within Dc which cannot 
be mapped to each other by this transformation. 

In section |2j I review the Ginsparg- Wilson chiral symmetry, to give the no- 
tation and tools which I will use in section [3] to construct the new Ginsparg 
Wilson operators and demonstrate that, with sufficiently smooth gauge fields, 
they are exponentially local. I test the locality numerically for a few examples 
in section H] and, after outlining some generalisations to this work in section 
[5j I conclude in section [6j 



2 Ginsparg- Wilson chiral symmetry 

2. 1 Introduction 

In this section, I review the demonstration that the Ginsparg- Wilson equation 
implies an exact lattice chiral symmetry and gives a lattice topological charge 
satisfying an index theorem. This section is based on the work of Martin 
Liischer [I]. 

Consider the following "chiral" transformation of the fermion fields: 
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which for small enough a can be written as 

if)'. = il) i + «7 5 (s - ^aRD^J ipi 

i J ' i =^i + (s - ^aDR^j 75. (3) 
If there is a mass-less fermionic action 

N f 

S f = J2AD^ (4) 
i=i 

where D is 75-Hermitian (D' = 75D75), then it is trivial to demonstrate 
that the action is conserved under the chiral transformation if and only if the 
general Ginsparg Wilson relation is fulfilled 

£75S + S l5 D = a^D(R l5 + r y 5 R)D. (5) 

We can write R = Ra + Rc, where Ra = |(-R— 75-R75) and Rc = \(R+^R~i5) 
so that {Ra,^} = and [i?c, 75] = 0. Since Ra does not contribute to the 
Ginsparg- Wilson relation, we can simplify without any loss of generalisation 
by restricting R to those functions which commute with 75. Additionally I 
restrict S to 75-Hermitian operators. By taking the Hermitian conjugate of 
equation (EI), it is clear that, with these restrictions, R must be Hermitian. It 
is now possible to simplify equation ([5]) so that it reads 

D*S + = aD^RD. (6) 



This can be associated with the continuum chiral symmetry because in the 
continuum limit (for suitable S and R, i.e. they are local and with the correct 
form in the continuum) it reduces to 75-D + D75 = 0. The Ginsparg- Wilson 
relation is usually defined with S = 1, and S can easily be absorbed into the 
definition of D and R. Of course, the Dirac operator must also satisfy the 
usual criteria for a suitable Dirac operator, namely that it is exponentially 
local (according to the definition D xy < ae~^ x ~ y ^ for positive a and /3), that 
the Dirac operator has no fermion doublers, and that expanding the Fourier 
transformed operators in the lattice spacing, a, gives 

D(p) =m 7At (^ + 4r b ) + 0(a 2 ) 
S(p) =1 + 0(a) 

R(p) =0(1), (7) 

where T b are the (Hermitian) generators of the gauge group, and = A h ^T h 
represents the gauge fields. It can be shown that a Ginsparg- Wilson Dirac op- 
erator correctly resolves the U(l) anomaly by noting that the fermion measure 
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is not invariant under this transformation: 



dip' dip' = det 



3 «75 



[S-^aRD) 



det 



e a 



S-^aDRy 



dip dip, 



which leads to a definition of a lattice topological index [4120] 

Q f = Tr( l5 S-^a l5 RD). (9) 

As an aside, I note that, in the continuum limit, the Yang Mills action can be 
written as [22] 

-F^ = c s Tr (S - -aRD) + constant, (10) 

and the electromagnetic field tensor as [23] 

F lu , = c f Tr a^(S-^aRD), (11) 

where c s and Cf are normalisation constants. This means that all the ele- 
ments of QCD and the electro- weak Lagrangians can be constructed from any 
Ginsparg- Wilson operator. 

R can, of course, be trivially constructed algebraically for any lattice Dirac 
operator (assuming that its inverse exists), 

R=s h + ^- (12) 

but only a few possible Dirac operators will give the local R needed for the 
lattice chiral symmetry. It is trivial to show that inserting the Wilson Dirac 
operator gives a non-local R since the inverse of the Wilson operator is non- 
local and has additive mass renormalisation 



One exact solution to the Ginsparg- Wilson equation, with local R, has been 
known for over ten years. The mass-less overlap operator [22], given by 



D = l+ T6 e(75D w ), (13) 

where e is the matrix sign function, satisfies the Ginsparg- Wilson equation 
with R = 1. D\y can in principle be any valid lattice Dirac operator with 
the correct continuum limit, no fermion doublers, and a suitably chosen mass 
term (e.g. in equation ffl4l) . any m satisfying m c < m < 2, where m c ~ is 
the Wilson critical mass, will suffice, although in practice m should be tuned 
to improve the locality of the operator). For the purposes of this work, I will 
use the simplest possibility, the Wilson Dirac operator, which I shall write as 



Dw 



~£M^ + ^)-^]-m. (14) 
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<9 M and d* are the forward and backward lattice Dirac operators, defined as 

d^{x)) = e iaA ^ {x+a ^^{x + ajj) - ijj{x) 
dl(iP(x)) = *P{x) - e- iaA ^ x - a ^^{x - a/2). (15) 

Thus, taking the Fourier transform to obtain the momentum representation 
of D w , I obtain 



D w (p) =«Tm^ + W ~ 



m 



ff^ — _ 1 e i(p,x) ^iaA^x+a^) e iap^ _ e ~iaA^{x-a^) ^iap^ 

2 ~ 



\Y = — — e^'^e^'^ f e iaA v(. x + a 2 ) e ia Pn _|_ e - iaA n(y- a % ) e ~ ia P^ 
Note that in the free theory, 

Hp = sm(ap^) 



Using equation (j!6p . the momentum representation of the overlap operator 
can be expressed as 



B(p) 

B(p) =(m - Wf + HJIp + {(W - m)H fJi - H^W - m)} . (18) 

It is clear that, because the term inside the square root is real and greater 
than or equal to zero, for gauge fields which are sufficiently smooth and where 
DyyDw does not have an exact zero eigenvalue, the Fourier representation of 
the overlap operator is an analytic function of the momentum, and that the 
momentum is bound — 7r < pa < ir. From the Paley-Weiner theorem [26J , this 
is enough to ensure exponential locality. Numerical experience has shown that 
the overlap operator remains exponentially local even when DyyDw has a zero 
eigenvalue [27]. Dq{p) is only zero at p = 0, meaning that it has no unwanted 
doublers. It is also clear that it has the correct continuum limit. The overlap 
topological charge is 

Qo = -^Tr (e( 75 iV)) (19) 
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2.2 The eigenvalue spectrum of generalised Ginsparg- Wilson operators 

I now demonstrate that those chiral Dirac operators satisfying both [D, D*] = 

and the Ginsparg- Wilson equation have an eigenvalue spectrum on a curve 
in the complex plane. This condition, along with 75-Hermiticity, guarantees 
that the eigenvalue spectrum is symmetric under reflection in the real axis. 

1 start with the Ginsparg Wilson equation 

D + D ] = D ] RD. (20) 

Using a spectral decomposition, which is valid if D and share eigenvalues, 
i.e. [D, £>t] = 0, I obtain 

2Re(A) = |A| 2 (/0, Rip). (21) 

It is clear that equation (1201) implies that 

[D,D*](1 -RD) =D^[D,R]D. (22) 

Thus, if [D, D'] = then [D, R] — and ip is also an eigenvector of R. This 
enables me to write 

2Re(A) = |A| 2 i?'(A), (23) 

where R' is some real function of the eigenvalue. Thus, in this case, the eigen- 
value spectrum is constrained to a curve in the complex plane. 

In the general case, the eigenvalues of a Ginsparg- Wilson operator will not lie 
on a curve. For example, we can consider a Dirac operator D, defined by [28] 

T^DT = D GW , (24) 

where D GW e D G is some Ginsparg- Wilson operator obeying the relation 

D G wl5 + l5D GW = D GW j 5 R GW D GW , (25) 

and T is some local operator which commutes with 75, satisfies the correct 
continuum limit and whose inverse both exists and is local. In the context 
of the renormalisation group and the fixed point action, T is equivalent to a 
change in the blocking procedure used to modify the lattice spacing. Given 
that T is invertible, D will satisfy a Ginsparg- Wilson equation 

D l5 + l5 D = D l5 TR GW T^D. (26) 

For every vector satisfying D G w \ip) = there will be a vector T~ l which 
is a zero eigenvector of D. Otherwise there will (in general) be no relation 
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between the eigenvalues and eigenvectors of D and Dqw unless T commutes 
with Dqw or is unitary. Every particular Ginsparg- Wilson operator is related 
to at least one operator within the class &c by some transformation T. The 
index of the Dirac operator is unchanged: 

Qt = -^Tr RD = -^Tr TRqwT^ ~ x Dqw CO _1 = Qgw- (27) 

Not every operator within 2)c can be mapped to every other operator in Dq- 
To move from one Dirac operator to another which commutes with it, it is 
necessary to use a transformation T which commutes with both Dirac opera- 
tors, and both R functions. In some cases R might have some zero eigenvalues, 
and the transformation needed to increase the number of zero eigenvalues of 
R is not invertible. 



3 Additional solutions to the Ginsparg- Wilson equation 



Consider the Dirac operator 



D r 



l + r 



756(75.0^), 



(28) 



where r[x] is some real, positive, and analytic function (which, as I shall 
demonstrate, is enough to ensure locality), and which satisfies 



r[±l] = 1. 



(29) 



The condition given in equation ff29|) ensures that, for a suitable choice of Dw, 
such as the one given in equation (TH1) . the operator has the correct continuum 
limit and no doublers. It is trivial to show that because r commutes with both 
75 and e, D r is 75-Hermitian. Furthermore, it satisfies the Ginsparg- Wilson 
relation with 



S =1 



R r 



D 



(2 + r 



^ (e7s + 7se)] (e7s + 7s<0) 
2 + (e7 5 + 75e)^ l( e T5 + 75e 



l + r 



|(e75 + 7se) 



|(e75 + 75e) (e7s + 75e 



(30) 
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From equation (flBI) . I obtain 



D r (p) =1 + r [x] (7 M ff M + W - m) ' 



1 

X= 2 



£(p) 

(7m#» + W(p) - m)-^— + -=— (- 7 „#» + W(p) - m) 

VB(p) v B(p) 

(31) 



It is clear that if the gauge field A and the function r are analytic and DyyDy/ 
has no zero eigenvalues, then the Fourier representation of the Dirac operator 
is analytic. Therefore, using the same argument as for the overlap operator, 
it is exponentially local. It is also clear, given equation f[2"§j) . that these Dirac 
operators have the correct continuum limit and that the doublers have infinite 
mass. In the momentum representation, R is given by 

5 2 + 2r [x] x 

AV) = l + (r [x]) 2 + 2r [x]x' (32) 

There is potentially a pole in R r at 1 + r[x] 2 + 2r[x]x = or, equivalently, at 
r = —x + \Jx 2 — 1. Since r and x are both real, — 1 < x < 1 and r is positive, 
the pole cannot contribute for any i / -1. Expanding R r around x = — 1 
gives 



Rr 



2(1 - r'[-l])(l + x) + |(1 + x) 2 (2r'[-l] - r"[-l]) + 0((1 + x) 3 ) 
2(1 + x) + |(1 + x) 2 (2r'[-l] 2 + 2r'[-l] + r"[-l]) + 0((1 + x) 3 ) ' 

(33) 



Therefore, R r has a smooth limit to x = —1 and does not diverge. B(p) is 
always real and positive (again assuming analytic r, A and that D^Dy/ has 
no exact zero eigenvalues), so R r is analytic, and R r will (at worst) fall of 
exponentially with distance. Thus both the Dirac operator and the Ginsparg- 
Wilson function R r are exponentially local, and this operator should be a 
suitable lattice Dirac operator with chiral symmetry. Given my experience 
with the overlap operator, I do not expect locality to break down when Dy^Dw 
has an exact zero eigenvalue. 

I now need to demonstrate that these operators have a well defined topological 
charge; which is most easily done by comparing the eigenvalue spectrum of 
this operator with that of the overlap operator. I write the non-zero modes of 
the Hermitian overlap operator as \ip+) and 1^-)) where the zero modes are 
|^o) and the unpaired eigenvectors of the overlap operator with eigenvalue ±2 
are Since D Q and D r commute, it is easy to show, that these eigenvectors 



9 



are also eigenvectors of 7s-D r and R r . Additionally, it can be shown that 



(i/; + \D r \^ + )=-(^\D r \^ 

R r |-0O> =Rr 1^2) = 1 
# r |V+> =i2r|V->- 



(34) 



Hence, 



--Tr 75-Rr-Dr 




75 A? 



(35) 



Therefore these Dirac operators will have the same zero modes and topological 
charge as the overlap operator. 

I note that writing the operator as r(cos8)e ld , which I have done in equation 
(1281) . is only possible because [Do, D Q ] = 0. It would not be possible to perform 
a similar decomposition with, for example, the Wilson operator to shift the 
Wilson eigenvalues to a closed curve on the complex plane. 



4 Numerical tests 

To test the locality of the Dirac operator D r and the Ginsparg- Wilson function 
R, I measured 



where <f)(x) is a point source. I then plot the mean value of the locality L(\x — 
y\ = d) against d, the distance between x and y in lattice units. 

The operator with r = l/\x\ mimics the continuum operator in that the eigen- 
value spectrum is on two lines — one along the imaginary axis containing the 
physical modes, and one at a constant mass (in the continuum limit an infinite 
mass), containing the doublers. However, because this operator is impracti- 
calj_j I shall use an approximation to this function. I therefore set r as the 
Chebyshev approximations (over the range 0.1 < x < 1) of order n — 1 of 

1 Since and W do not commute in the non-Abelian theory, I have only been to 
write the operator in a simplified form in the free or Abelian theory. To calculate 
the operator exactly would require calculating eigenvalues of the overlap operator 
to a high precision and using some rational or polynomial approximation to \x\ to 
simulate the rest of the eigenvalue spectrum, which is possible but impractical. It 
is also not clear if this operator is local: r is not analytic; however we can construct 



L D (x,y)=(<l>(x)\D r \<f>(y)), 
L R (x,y)=(<j ) (x)\R r \<j ) (y)) 



(36) 
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Fig. 1. The eigenvalue distributions of the various Dirac operators used in this study. 

the function / = l/y/x^, for n — 1, 5, 9 and 13, normalised so that it satisfies 
equation fl29|) . This gives me four different Dirac operators, Di, D§, Dg and 
-D13. The (theoretical) eigenvalue distributions for the operators D n are shown 
in figure [TJ Di is the overlap operator. The kernel of the matrix sign function, 
D w , is as defined by equation (TH1) . but with two levels of stout smearing [3U] 
at parameter 0.1, and with m = 1.5. I tested the locality of these operators on 
configurations from a 12 3 48 Dynamical overlap ensemble, with lattice spacing 
a ~ 0.13fm. I apply the Dirac operator D n or Ginsparg- Wilson function R n 
to a unit source at x, and calculate its projection onto another unit source 
vector at y, averaging over x, y and configurations. The results are shown in 
figures [2] and [3j 

It can be seen that in all cases both the Dirac operator and the Ginsparg- 
Wilson function are exponentially local, with the same rate of decay as the 
overlap operator. They are all constrained so that L D < ct^e - ^ - ^ and L R < 
a^e _/3 ' 1_?/ '. The value of (3 seems to be independent of n, and (we might 



successively better approximations to r all of which will be local. I also note that 
because R\/\ x \ = for the physical modes (it is not zero for the eigenvalues cor- 
responding to the doublers; so that the Nielson-Ninomiya theorem is still avoided) 
the left and right chiral projectors are identical for the physical modes. Thus if this 
operator is local, it might be possible to use it to circumvent the CP and T viola- 
tions described in [29]. However, a full discussion of this interesting topic is beyond 
the scope of this paper, and deserves a full treatment in a future work. 
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conjecture) seems to be general across all possible functions r. The value of a 
increases with increasing n for the Dirac operator (although there is very small 
change for the Ginsparg- Wilson function), and does depend on r. It is possible 
to suspect that other forms of the Ginsparg- Wilson Dirac operator could have 
improved locality compared to the overlap operator[f] It can be concluded 
that these Dirac operators and Ginsparg Wilson functions are, as expected, 
exponentially local. The various small wiggles in the curves are caused by the 
geometry of the lattice and the breaking of rotational symmetry. 



5 Generalisation of the Ginsparg- Wilson Dirac operator 



The proposed operator given in equation (|28|) can be easily generalised further. 
One possibility is to multiply by a function q: 



D 



rq 



(7se + e7s) 



1 + r 



(7 5 e + e7 5 ) 



75e 



(37) 



where we again have the constraints that q[x] must be positive, analytic, and 
q[l] = q[—l] (for the sake of being definite, I shall use g[±l] = 1, although the 
value can easily be absorbed into the fermion renormalisation constant), and, 
to give the correct continuum limit, restricted to q — 1 + 0(a). This leads to 
a Ginsparg- Wilson function 



R 



rq 



q[x\ 



7 R r \x\ 



(3* 



and it is clear that, with q as specified, this contains no poles in the Fourier 
representation, so R rq will be local. We can, of course, generalise further, and 
consider an operator 



D 



rqth 



t 



^(75e + e75) 



h q 



2(7se + e7 5 ) 



1 + r 



2(75e + e7s) 



75 e 



(39) 



where h(z) is analytic, positive except for z = and satisfies h(0) = 0, h(2) > 
and h(zY = h(z^); and t obeys the same conditions as q and r. Finally, I can 
choose 5* 7^ 1 in equation (Q, which will modify the above Dirac operator to 
give: 



D 



rqthS 



S( T5 e)1 



2-(7se + e75) 



h q 



2-(7se + e75) 



1 + r 



2"(75e + £75) 



75e 
(40) 



2 For the same kernel. It is known that using an improved kernel can improve the 
locality [31], but this is a separate issue. 
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where S(z) is an analytic function satisfying S(—l) = S(l) = 1 and S(z^) = 
S(zY . Once again, from the Fourier representation, it can be shown that D rqt hs 
and Rrqths are both exponentially local and have the correct continuum limit. 
Fujikawa's operators [TTJ are members of this class, with r = q = S = t = l, 
h(z) = z l ^ 2k+l \ and a particular form for the the kernel operator Dw which 
gives the correct continuum limit. 

Kerler [13] considered, at first, a alternative form of lattice chiral symmetry 

lb D + D lb V = 0, (41) 

where V is unitary and 75-Hermitian. This relates to the standard Ginsparg- 
Wilson relation with the substitution of variables V = 1 — RD = —j^D. 
He later extended equation (14"T|) to include some possible S 7^ 1 forms of the 
Ginsparg- Wilson equation (following the transformations outlined in [29|), 
and, by considering Dirac operators which are functions of V and using a 
spectral decomposition, he showed that operators of the form 



D K = -i(G(V)G(V)) 1/2 H i^:(V- 1/2 - V 1/2 )W 



(42) 



where H(—z) = —H[z\ (with the expansion around z = giving the correct 
continuum limit; for example using Dw as the kernel of the matrix sign sign 
function, H[z] = z + 0(z 2 )), W(—l) 7^ 0, both functions are Hermitian and 

G(V) =((l-8 k ) + 8k V)/N 
G =(s k + (1 - s k )V)/N 



N= x jl- 2s k {\ - s k ){\ - 1 -{V + Vt), (43) 

satisfy a lattice chiral symmetry. It is clear that Kerler's solution (with V = 
756) and my own have a certain similarity, although they were derived from 
different approaches: his from attempting an algebraic solution of the eigen- 
value equivalent of the Ginsparg- Wilson equation; mine from desiring to test 
a seemingly general property of chiral Dirac operators. His solutions are a par- 
ticular form of my most general D rqt hs solution, specifically with r = 1 and 
with a particular form of S and t. Kerler did not discuss in detail the locality 
of the Ginsparg- Wilson function, which is crucial in determining if we truly 
have a lattice chiral symmetry. It is unclear whether these generalisations offer 
any advantage over the operator given in the previous sections!"^! 

Finally, it is worth spending a moment considering whether it is necessary to 



3 One possibility might be to absorb a perturbative expansion of the fermion renor- 
malisation into the Dirac operator; or to otherwise reduce the higher order lattice 
artifacts. 
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use 756 to construct these operators, and not some other unitary operator u. 
There are several ways in which we can construct u: 

(1) u = C/VC^C; 

(2) u = -^-C for C not ( ant i-) Hermit ian; 

(3) u = e A , for ant i- Hermit ian A. 

u, of course, would have to be local and have the correct continuum limit. 
I have already discussed the first option. For the second option, forcing u 
to be local, having the correct continuum limit and being 75-Hermitian will 
place considerable restrictions on the possible choices of C (one option is 
C — 1 + Dqw, where Dew is another Ginsparg- Wilson operator; but if we 
choose to use a Ginsparg- Wilson operator constructed the matrix sign function 
this will revert to another form of equation (|40p . For the third option, it is 
not obvious how to construct suitable operators which are both 75-Hermitian 
and free of doublers; for example 275(1 — e Mr// ' /3 ' 75Dw ) is not 75-Hermitian, 
and 1 — e D ~ l5C>l5 has doublers. Thus using the matrix sign function remains 
the only currently known possibility to construct local chiral lattice Dirac 
operators free from doublers. 

One can also, of course, transform the Dirac operators using the method out- 
lined in section [2721 By choosing T and T* which commute with D, we can map 
the eigenvalue spectrum to any closed curve (which passes through the origin 
to get the correct continuum limit and through two to remove the doublers). 
The only restriction with the mapping is caused by the number of zero modes 
of R; but since R rq ths has an arbitrary number of zero modes depending on 
the choice of functions, we can map to any possible R and hence D. Thus 
the eigenvalue spectrum of any Dirac operator within &c can be mapped to 
the eigenvalue spectrum of at least one of the operators D rqthS by a suitable 
transformation. The eigenvectors of D^^s are determined by Dw, and can be 
freely modified (within certain constraints) by adjusting Dw- Given that the 
Dirac operators can be defined in terms of their eigenvalues and eigenvectors 
spectrum, it is plausible, if not yet proven, that all Dirac operators within ®c 
are Ginsparg- Wilson operators. 



6 Conclusion 

I have demonstrated that the overlap operator is just the simplest member of 
a class of chiral Dirac operators by constructing additional Ginsparg- Wilson 
operators. I propose that any lattice Dirac operator which has 

(1) The correct continuum limit with no doublers; 

(2) Eigenvalues which lie on a closed loop in the complex plane, symmetric 
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under reflection in the imaginary axis, and single valued with respect to 
the angle from the center of the Ginsparg- Wilson circle; 

will be exponentially local (in a sufficiently large volume and sufficiently 
smooth gauge fields), will satisfy a lattice chiral symmetry with exponentially 
local (or better) Ginsparg- Wilson function R with an exact index theorem, 
will be 75-Hermitian, and will thus be a suitable lattice Dirac operator. Of 
course, these additional Dirac operators are more expensive to simulate than 
the overlap operator while it is not clear that they have any benefits over 
the overlap operator. Thus it is unlikely that they will have any more than 
theoretical interest. 

Furthermore, I have shown that chiral Dirac operators fall into certain homo- 
topy classes determined by the number of zero modes of the Ginsparg- Wilson 
function R. In the continuum limit, all these Dirac operators will reduce to 
the same, universal, operator, and it is most unlikely that this division has any 
physical relevance. However, it seems likely that the perfect action (the lattice 
action with no discretization errors) will have to fall into the same homotopy 
class as the continuum Dirac operator. 
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